Charge-Hall effect driven by spin force: reciprocal of the spin-Hall effect 
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A new kind of charge-Hall effect is shown. Unlike in the usual Hall effect, the driving force in the 
longitudinal direction is a spin force, which may originate from the gradient of a Zeeman field or a 
spin-dependent chemical potential. The transverse force is provided by a Berry curvature in a mixed 
position-momentum space. We can establish an Onsager relation between this effect and the spin- 
Hall effect provided the spin current in the latter is modified by a torque dipole contribution. This 
remarkable relation leads to new ways for experimental detection of spin accumulation predicted by 
the spin Hall effect. 
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The generation of an electric current in the transverse 
direction of an electric field is known as a Hall effect. The 
transverse force is usually provided by the Lorentz force 
from a magnetic field, but can also arise from spin-orbit 
interactions in magnetic materials |l| or Magnus forces on 
vortices in superconductors |2J. Apart from such ordinary 
and anomalous Hall effects, other variants, such as the 
thermal Hall effect 3], can manifest when the longitudi- 
nal electric field is replaced by other type of electromotive 
forces like a temperature gradient. Likewise, a longitu- 
dinal electric field can also cause non-electrical response 
in the transverse direction, such as in the spin Hall effect 
where a transverse spin current may be generatedka,|5j,|6j. 

Recently, extensive theoretical workJ3,|3,y,|lfJ has es- 
tablished the importance of intrinsic origin of the anoma- 
lous Hall effect in ferromagnets. Unlike the traditional 
mechanisms of skew scattering [ll|, Berry phase in the 
Bloch bands due to spin-orbit coupling can lead to a 
Hall current depending only on the equilibrium part of 
the carrier distribution function, fn a similar fashion, 
intrinsic spin-Hall effect has been proposed for paramag- 



netic materials 
est recently [l4l 




1 131 . which ha s ge nerated much inter- 
El [H El H Hail miH in associ- 
ation with electrical manipulation of spins for spintronics 
applications 24, 25, 26] . 

In this Letter, we propose a novel charge Hall effect 
driven by a spin force along the longitudinal direction. 
This spin force can be provided by the gradient of a Zee- 
man field or a spin dependent chemical potential. In the 
intrinsic regime considered here, the transverse force is 
provided by a Berry curvature in the mixed position- 
momentum space, also stemming from spin-orbit cou- 
pling in the band structure. As shown below, there is not 
a longitudinal electric current accompanying the spin- 
force in this intrinsic regime, so the spin force should not 
be interpreted as a variant electromotive force. There- 
fore, this effect is essentially different from the conven- 
tional thermal Hall effect, which is known to be just a 
usual Hall effect driven by a variant electromotive force. 

Another goal of this Letter is to construct a "recipro- 
cal" version of the intrinsic spin-Hall effect [12lll3lll4ll5L 
HElElElliillillilEl which is at present being 



in active investigation. Since a longitudinal charge force 
(electric field) may drive a spin-Hall current through the 
spin-orbit interaction, one naturally expects that a spin 
force can also induce a charge-Hall current. To our sur- 
prise, we discovered that an exact Onsager relation be- 
tween these intrinsic Hall effects cannot readily be estab- 
lished; only when the spin current is modified by includ- 
ing a torque dipole term, is the spin-Hall conductivity 
in response to an electric field found to correspond to 
our charge-Hall conductivity in response to a spin force. 
Interestingly, it is this modified spin current that is re- 
sponsible for spin accumulation at a sample boundary as 
shown recently based on the spin continuity equation in 
a semiclassical theory |l4J. Our finding on the Onsager 
relation further points to the importance of this modi- 
fied spin current. Spin accumulation due to spin Hall 
effect can thus be tested through the Onsager relation 
by a measurement of our charge Hall effect in response 
to a spin force. 

To construct the theory, we consider the conduction 
electrons in semiconductor quantum wells. In this two- 
dimensional electron gas (2DEG) system, spin-orbit in- 
teraction arises from the quantum well asymmetry in the 
growth (z) direction and has a standard Rashba form|27l] 
V so = (a/h)(a x k) ■ z, where a is the Pauli matrix vec- 
tor and k is the 2D wave vector in the x-y plane. The 
Rashba coefficient, a, can be tuned over a wide range by 
a vertical electric field. Taking into account the kinetic 
energy and a Zeeman field normal to the heterostructure, 
we can write the k • p Hamiltonian as 



H = jk + V so - \io z 



(1) 



where 7 = H 2 /2m* with m* being the band effective 
mass, and /i is the Zeeman field (in unit of energy), which 
is assumed to be weak in this paper. The orbital effect of 
the magnetic field can be taken into account in our semi- 
classical formalism, but will be neglected in this work for 
simplicity. To introduce a spin force, we consider a non- 
uniform Zeeman field as in the Stern-Gerlach experiment 
for separating spin-up and spin-down electrons. We thus 
assume /i = /io + fiix, where /jg describes the average 
strength of the Zeeman field while F s = (i)d/j,/dx = \[i\ 



gives the spin force in the x (longitudinal) direction. The 
Zeeman field can be created by exchange interaction with 
the moments of doped magnetic ions, which can be po- 
larized by a weak magnetic field (with negligible Lorentz 
force) at low temperatures. An inhomogeneous Zeeman 
field can be produced either by a magnetic field gradient 
or a temperature gradient. Alternatively, one may con- 
sider using the spin-dependent chemical potential gradi- 
ent present near an interface with a ferromagnetic mate- 
rial. 

To access transport properties, we adopt the formal- 
ism of semiclassical wavepacket dynamics 28], which is a 
powerful tool for studying the influence of slowly varying 
perturbations such as the spin force term in Eq.(l) on 
the dynamics of Bloch electrons. Consider a wave packet 
centered at r c and with a narrow distribution around the 
mean wave vector k. Then the semiclassical equations of 
motion for r r and k are |28| 



17 = 



1 de n 
h dk ' 



l<9Ae„ 

h 9k 

l dAe„ 

fi dr r 



n h 



n: 



r c - ^r • k, 



a rk ■ k, 



(2a) 
(2b) 



where n is the band index, and e n is the local band en- 
ergy obtained by diagonalizing the Hamiltonian (1) with 
the position operator r identified with the wavepacket 
center r c . The result is £i,2(k, r c ) = jk 2 =p y/(i% + a 2 k 2 
with fi c = fj,o + [i\x c . In the second term in Eqs.(2), 
Ae„ is the energy correction in the gradient of the Zee- 
man field. For the Hamiltonian (1) we get Aei(k, r c ) = 
Ae 2 (k, r c ) = a 2 k y ^i/2(^ 2 + a 2 k 2 ). Finally, the last two 
terms in Eqs.(2) denote Berry-curvature corrections to 
the orbital motion of the electron. They are defined by 



K r ) afJ ^rt 



c,3 _ 



/ dun dun _ , du n . du n . 
dk a dxp dxj3 dk a 



, (3) 



where u n is the eigenstate of the nth band. The other 
Berry curvature tensors are defined similarly. One can 
see that the equations of motion involve Berry curvatures 
between every pair of parameters in combined real and 
momentum spaces. Recent studiesJ3,Hla| of the anoma- 
lous Hall effect in ferromagnetic materials have revealed 
the fundamental role of the Berry curvature J7 kk in the 
k space. As we will see below, it is the Berry curvature 
in the mixed position-momentum space, i7 kr , which has 
not attracted attention in previous work, that determines 
the charge-Hall effect driven by a spin force. 

In most cases, calculation of the Berry curvatures is an 
involved work. For real materials with complicated crys- 
tal and ferromagnetic structures, first-principle calcula- 
tions are needed. Fortunately, in the present two-band 
case, the Berry curvatures can be obtained analytically. 
For the bottom band with energy ei, the non-zero com- 



ponents of the Berry curvatures are 

o?k y [i 



n 



fi'jy ■ "■■ ii 



2A 3 . 



fi**'* 



2A? 



i aky,x 



a 2 k x ^i 
2A? 



^ ' (4) 

where A c = y/fj. 2 . + a 2 k 2 . The Berry curvatures for the 
top band are different from the above by a sign. Consid- 
ering these non-zero Berry curvatures, and keeping terms 
up to first order in spin force, we get k x = —de n /hdx c 
and k y = 0. The orbital motion in the nth band is given 
by 



(!-(#■>» + 



dAe n 
hdk x ' 



Vc 



dAe 71 
hdky 



fL ky ' x v* + n^' k " 



de n 
hdx c 



(5a) 
(5b) 



where v 1 ^ = (l/H)de n /dk a is band group velocity. The 
longitudinal and transverse currents are obtained by av- 
eraging x c and y c with the distribution function, re- 
spectively. The full distribution function consists of an 
equilibrium part and a noncquilibrium part (maintained 
by a competition between fc^-drift and scattering relax- 
ation). As for the anomalous Hall effect and spin Hall 
effect, there is an intrinsic contribution from the equilib- 
rium part of the distribution. The extrinsic contribution 
from the non-equilibrium part will be small if scatterings 
are strong enough to maintain the distribution near to 
equilibrium and weak enough (compared to the Rashba 
splitting in energy scale) so that interband couplings are 
small. In this work we will focus attention to the intrin- 
sic regime in order to make contact with similar work for 
the spin Hall effect. Using the equilibrium distribution, 
we find that the longitudinal current vanishes. However, 
the intrinsic Hall current in the y direction is not zero, 
and is given by, using Eq.(5b), 



Jy-^ 



-R 



d 2 k 



f D. k v x v x 



(2n) 2 

d 2 k eha 2 k, 
(2tt) 2 4A 3 



(6) 



ihvl ~ f2V X 2 ) 



rr cs F 



where /„ is the equilibrium Fermi-Dirac distribution 
function, and A = (uQ + o^fc 2 ) 1 / 2 . We note that the con- 
tributions from the second and last terms in Eq. (5b) can- 
cel each other under the k integral, thus only the Berry 
curvature U ky ' x contributes to the charge-Hall current. 

Equation (6) is a central result in this paper. Figure 
1 summarizes the behavior of the charge-Hall conductiv- 
ity. Panel (a) is a (colored) contour plot of this quantity 
in the parameter plane of electron density and Rashba 
coefficient. The (blue) region below the dashed line cor- 
responds to the case of single band occupation, where the 
charge-Hall conductivity is small. At high electron den- 
sities, cjy X remains large within a wide range of a. Panel 
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FIG. 1: Charge-Hall conductivity (in unit of e) in GaAs quan- 
tum well system as a function of (a) Rashba coefficient a 
and the electron density iV e (contour plot) for T = and 
/io = O.OlmeV, (b) Rashba coefficient for different values of 
Ho, and (c) temperature. 



(b) shows the dependence on the Rashba coefficient for 
a set of Zeeman energies. In these panels, the tempera- 
ture is set to zero. Panel (c) shows that the charge-Hall 
effect is robust against thermal broadening, a property 
stemming from the intrinsic nature of crfj*. For estimate, 
we consider a GaAs quantum well sample with electron 
density iV e ~ 5 x 10 n cm~ 2 and a typical Rashba coeffi- 
cient of a ~ 20meV-nm. A spin force corresponding to 
fii ~ lmeV//xm (one millivolt of Zeeman energy change 
over a micron) can yield a charge-Hall current on the 
order of ImA/cm, implying that the predicted effect is 
observable. 

Now we turn our attention to a basic and enlighten 
question: is there an Onsager relation between our 
spin force-driven charge-Hall effect and the electric field- 
driven spin-Hall effect? In the presence of an electric field 
along the y direction, the Hamiltonian is now 



H = 7 fc 2 



V m 



A*oc 2 



eEy, 



(7) 



where we have removed the spin force but kept the uni- 
form part of the Zeeman field. Based on the spin current 
operator, defined as the symmetric product of the veloc- 
ity and spin operators, h(s z x + is z ), one can obtain a 
spin current in the transverse (x) direction either from 
the Kubo formula or the semiclassical wavepacket for- 
malism. In the intrinsic regime, this spin current for the 
Rashaba model is given by 






-E 



d 2 k eh 2 a 2 k 2 . „ „ . 
#(/i-/ 2 ). 



(2tt) 2 4m*A 3 



(8) 



At zero temperature and up = 0, this yields a universal 
spin Hall conductivity EH EE E3, EE E3, Q E3l , <% = 



—(Ty% = —e/8ir, which does not depend on the Rashba 
coefficient. Unfortunately, this universality is absent in 
the spin force-driven charge-Hall conductivity a°| [see 
Eq.(6) and Fig.l]. In fact, the expression for o~™ in Eq.(6) 
can be reduced to the expression for <j x c y in Eq.(8) only 
when the band group velocity vf 2 = hk x /m*^a 2 /hA 2 is 
replaced by its first term hk x /m*; the neglected term will 
cause a finite difference for non-zero Rashba coefficient. 
Thus contrary to our intuitive expectation, a full Onsager 
relation between these two kinds of Hall conductivities 
fails to be established. 

To understand the origin of this violation of the On- 
sager relation and to search for a remedy, we notice that 
the spin-current operator is not really an observable in 
conjugation with the spin force. To illustrate this point, 
let us put the charge force (the electric field) and the spin 
force in an equal footing by including both of them in the 
perturbative Hamiltonian 



SH 



Fid, -Fod 



2"2, 



(9) 



where F\ = E and Fi = F s are the generalized forces ap- 
plied on the charge and spin degrees of freedom, whereas 
d\ = —ey and g?2 = s z x are the corresponding displace- 
ment operators. The response currents are obtained as 
expectation values of the generalized velocity operators 
in the perturbed state 



ch 
d\ 



-ey, 



(10a) 
(10b) 



where symmetrization between operators in Eq.(lOb) is 
implied. It is clear now that the current conjugate to the 
spin force (Fp. ) not only involves the spin current given 
from s z x [EjLUj, but also a torque dipole term s z x. We 
therefore introduce a modified spin current defined by 



*"?/ 



d 2 k 

(27T) 2 



Jn\Sz'£ \ SzX)m 



(ii) 



where the bracket indicates quantum mechanical average 
over the wavepacket made out of the nth band. Due to 
the spin-orbit interaction, besides driving the wavepacket 
in the k-space, the electric field will also give rise to a 
nonadiabatic correction to the spin wavefunction. This 
correction has been taken into account when constructing 
our wavepacket for each band. 

In terms of the wave packet center coordinate x c and 
velocity x c , we can rewrite the above expression as |l4| 



^ = E 



d 2 k 

(27T) 



T In 



<-.>n+ 5 w: 



Xc d& z)1 



(12) 
where (p x z ) n = {(x — x c )s z ) n is the x-component of the 
spin dipole in the wave packet. To first order in the elec- 
tric field, the time derivatives in the second and third 



terms can be replaced by —eE(d/Hdk y ), and the quanti- 
ties under the time derivatives can be evaluated at zero 
field as (s z )i,2 = ±MoS/2A, and {p s x ) 1 2 = — fia 2 fc y /4A 2 . 
It is then clear from symmetry that the last term in 
the above equation vanishes after the k integral. The 
wave packet velocity has a first order field correction 
due to Berry phase, x c = v„ + (e/H)E£ln y , where 
v® = de n /Hdk x is the band group velocity. The Berry 
curvature tensor r?^ k has been given in Eq.(4). Similarly, 
the wave packet spin also has a first order field correc- 
tion, given by E(s z )' n with (s z )' 12 = Tea 2 k x / '4A 3 . After 
these considerations, the modified spin-Hall current can 
be written as 



n 

e d 
h dk v 



d 2 k 

(2tt) 2 
(*#)„]■ 



fnK(s z )' n + -St* k y{s z ) n (13) 



Remarkably, we find that the contributions from the 
Berry curvature term and the spin dipole term com- 
pletely cancel each other. Thus only the first term in 
Eq.(13) contributes to the modified spin current, yield- 
ing a spin Hall conductivity 



d 2 k eha 2 k x 
(2tt) 2 4A 3 



(/i«i ~ /2U2 



(14) 



A comparison between Eq.(14) and Eq.(6) gives the 
Onsager relation 



VyU^o) =Vx C y(Vo)- 



(15) 



It should be pointed out that these two kinds of Hall 
conductivities have time-reversal symmetry T, which is 
absent in the conventional Ohm's law. Since the spin 
force and the charge-Hall current are odd under T, thus 
they can be related via T-invariant a yx . On the other 
hand, the spin current and the electric field are even un- 
der T; the resultant spin-Hall conductivity a x c y is also 
T-invariant. Due to this time-reversal symmetry, the 
general Onsager relation cr a p(fi>o) = Gfia{— Mo) reduces 
to Eg. (15) in the present context. We have also shown 
explicitly[29j that Eq.(15) remains to be true for the 4- 
band Luttinger model, and we believe that this is a gen- 
erally valid relationship. 

This remarkable Onsager relation is another central re- 
sult in this paper. As a consequence, the spin Hall con- 
ductivity (for the modified spin current) should behave 
the same as we depicted in Fig. 1 for our charge Hall 
conductivity in response to a spin force. All previous 
discussions on spin Hall effect measurement are based on 
spin accumulation on a sample boundary due to a bulk 
spin current; Ref.|l4| showed that it is the modified spin 
current that is relevant to such a measurement. The On- 
sager relation derived here provides another method for 
testing this modified spin Hall effect. 



In conclusion, we have shown that a spin force can 
drive a novel charge-Hall current through the spin-orbit 
coupling. The predicted charge Hall effect has a remark- 
able Onsager relation with a spin-Hall effect driven by 
an electric field. This Onsager relation is only attainable 
when the spin current is corrected by a torque dipole 
density. 
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